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Abstract: The solutions of elliptic problems with a Dirac measure in right-hand side are not 
H l and therefore the convergence of the finite element solutions is suboptimal. Graded meshes are 
standard remedy to recover quasi-optimality, namely optimality up to a log-factor, for low order 
finite elements in L 2 -norm. Optimal (or quasi-optimal for the lowest order case) convergence has 
been shown in L 2 -norm on a subdomain which excludes the singularity. Here, on such subdomains, 
we show a quasi-optimal convergence for the TP-norm, s ^ 1, and an optimal convergence in H l - 
norm for the lowest order case, on a family of quasi-uniform meshes in dimension 2. The study of 
this problem is motivated by the use of the Dirac measure as a reduced model in physical problems, 
for which high accuracy of the finite element method at the singularity is not required. Our results 
are obtained using local Nitsche and Schatz-type error estimates, a weak version of Aubin-Nitsche 
duality lemma and a discrete inf-sup condition. These theoretical results are confirmed by numerical 
illustrations. 

Key words: Dirichlet problem, Dirac measure, Green function, finite element method, local 
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1 Introduction. 


This paper deals with the accuracy of the finite element method on elliptic problems with a singular 
right-hand side. More precisely, let us consider the Dirichlet problem 


( p s) 


-A us 
u 5 


j XO in 
0 on dfl, 


where H c l 2 is a bounded open C 00 domain or a square, and 5 Xo denotes the Dirac measure 
concentrated at a point xq e R such that dist(xo, dR) > 0. 

Problems of this type occur in many applications from different areas, like in the mathematical 
modeling of electromagnetic fields [17]. Dirac measures can also be found on the right-hand side 
of adjoint equations in optimal control of elliptic problems with state constraints [ ]. As further 
examples where such measures play an important role, we mention controllability for elliptic and 
parabolic equations [9, 10, 21] and parameter identification problems with pointwise measurements 

[23]. 

Our interest in (Ps) is motivated by the modeling of the movement of a thin structure in a 
viscous fluid, such as cilia involved in the muco-ciliary transport in the lung [15]. In the asymptotic 
of a zero diameter cilium with an infinite velocity, the cilium is modelled by a lineic Dirac of force 
in the source term. In order to make the computations easier, the lineic Dirac of force can be 
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approximated by a sum of punctual Dirac forces distributed along the cilium [20]. In this paper, 
we address a scalar version of this problem: problem (Ps)- 

In the regular case, namely the Laplace problem with a regular right-hand side, the finite element 
solution u h is well-defined and for u e H k+1 (fl), we have, for all 0 ^ s ^ 1, 

\\u - u h \\ s ^ Ch k+1 ~ s \\u\\ k +i, (1) 

where k is the degree of the method [ ] and h the mesh size. In dimension 1, the solution u$ of 
Problem (Ps) belongs to fP(f2), but it is not H 2 (Q). In this case, the numerical solution and 
the exact solution us can be computed explicitly. If .To matches with a node of the discretization, 
Ufi = us- Otherwise, this equality is true only on the complementary of the element which contains 
To, and the convergence orders are 1/2 and 3/2 respectively in IL 1 -norm and L 2 -norm. In dimension 
2, Problem (Ps) has no H 1 (Q)-solution, and so, although the finite element solution can be defined, 
the IL 1 (D)-error has no sense and the L 2 (0)-error estimates cannot be obtained by the Aubin- 
Nitsche method without modification. 

Let us review the literature about error estimates for problem (Ps), starting with discretizations 
on quasi-uniform meshes. Babuska [ ] showed a L 2 (0)-convergence of order /i 1-e , e > 0 , for a 
two-dimensional smooth domain. Scott proved in [ ] an a priori error estimates of order 2 — |, 

where the dimension d is 2 or 3. The same result has been proved by Casas [ ] for general Borel 
measures on the right-hand side. 

To the best of our knowledge, in order to improve the convergence order, Eriksson [ ] was 

the first who studied the influence of locally refined meshes near xq. Using results from [24], he 
proved a convergence of order k and k + 1 in the VU 1,1 (D)-norm and the L 1 (D)-norm respectively, 
for approximations with a P^-finite element method. Recently, by Apel and co-authors [2], a L 2 (D)- 
error estimate of order h 2 \ ln/i| 3,/2 has been proved in dimension 2, using graded meshes. Optimal 
convergence rates with graded meshes were also recovered by D’Angelo [12] using weighted Sobolev 
spaces. A posteriori error estimates in weighted spaces have been established by Agnelli and co¬ 
authors [ ]. 

These theoretical a priori results for finite elements using graded meshes increase the complexity 
of the meshing and the computational cost, even if the mesh is refined only locally, especially if there 
are several Dirac measures in the right-hand side. Eriksson [13] developed a numerical method to 
solve the problem and recover the optimal convergence rate: the numerical solution is searched in 
the form uq + Vf, where uq contains the singularity of the solution and Vh is the numerical solution 
of a smooth problem. This method has been developped in the case of the Stokes problem in [20]. 

However, in applications, the Dirac measure at To is often a model reduction approach, and a 
high accuracy at To of the finite element method is not necessary. Thus, it is interesting to study 
the error on a fixed subdomain which excludes the singularity. Recently, Koppl and Wohlmuth 
have shown in [ ] quasi-optimal convergence for low order in L 2 -norm for Lagrange finite elements 

and optimal convergence for higher order. In this paper, we show in dimension 2 a quasi-optimal 
convergence in LP-norm, s ^ 1, and an optimal convergence in the case of low order. The L 2 -error 
estimates established in [19] are not used and the proof is based on different arguments. These 
results imply that graded meshes are not required to recover optimality far from the singularity and 
that there are no pollution effects. 

The paper is organized as follows. Our main results are formulated in section 2 after the Nitsche 
and Schatz Theorem, which is an important tool for the proof presented in section 3. In section 4 
another argument is presented to obtain an optimal estimate in the particular case of the Pi-finite 
elements. Finally, we illustrate in section 5 our theoretical results by some numerical simulations. 
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2 Main results. 

In this section, we define all the notations used in this paper, formulate our main results and recall 
an important tool for the proof, the Nitsche and Schatz Theorem. 

2.1 Notations. 



Figure 1: Domains Do and Dj. 

For a domain D, we will denote by || • \\ s ,p,D (respectively | • \ s ,p,D ) the norm (respectively the semi¬ 
norm) of the Sobolev space W S,P (D), while || • \\ St £) (respectively | • \ s ,d) will stand for the norm 
(respectively the semi-norm) of the Sobolev space H S (D). 

For the numerical solution, let us introduce a family of quasi-uniform simplicial triangulations 
Th, °f D and an order k finite element space Vjf cz JZq(D). To ensure that the numerical solution 
is well-defined, the space Vt is assumed to contain only continuous functions. The finite element 
solution ug e Vjf of problem ( P $) is defined by 

f Vu h 5 -\7v h = v h (x 0 ), Vv h e v£. (2) 

Jo 

We will also evaluate the ZP-norm of the error on a subdomain of D which does not contain the 
singularity, for s ^ 2, and, whenever we do so, we will of course assume the finite elements H s - 
conforming. We fix two subdomains of D, named Do and Di, such that Do czcz Di cc D and 
xo f Di (see Figure 1). We consider a mesh which satisfies the following condition: 

Assumption 1. For some ho, we have for all 0 < h ^ ho (see Figure 1), 

D™ n D? = 0, where D™ = (J T, 

TeT h 

7m, #0 


and DJ is the complement o/Di in D. 


2.2 Regularity of the solution us. 

In this subsection, we focus on the singularity of the solution, which is the main difficulty in the 
study of this kind of problems. In dimension 2, problem ( P $) has a unique variational solution 
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us e W^ :P ({X) for all p e [1, 2[ (see for instance [ ]). Indeed, denoting by G the Green function, G is 
defined by 

G{x) = ~^ lo g(M)- 

This function G satisfies — AG = 5 o, so that G(- — xo) contains the singular part of u$. As it is done 
in [3], the solution us can be built by adding to G(- — xo) a corrector term uj e solution of 

the Laplace ( problem 

—A uj = 0 in Q. 

uj = —G(- — xo ) on dfl. 

Then, the solution is given by 


(3) 


us(x) = G(x - x 0 ) + uj(x) = -— log(|a; - x 0 |) + uj(x). 

Z7T 

It is easy to verify that us f Hq(Q). Actually, we can specify how the quantity |[u<s||i lP ,n goes 
to infinity when p goes to 2, with p < 2. According to the foregoing, if we write us = G + uj, since 
uj e B' 1 (Q), estimating || 1 | i,p,r 2 when p converges to 2 from below (which will be denoted by p f 2 ) 
is reduced to estimate ||G||i iP) b, where B = B( 0,1): G e L P (I2) for all 1 =% p < oo, and using polar 
coordinates, we get, for p < 2, 


| G\ p lpB = £ |VG(x)| p dx = £ J£ (j-±y rdedr = (2 7 r) 1 ^£ r^dr 


(2tt) 


i -v 


2 -p 


Finally, when p / 2, 


11^.5 ||i, p, n 


V2tt V 2 - V 


(4) 


2.3 Nitsche and Schatz Theorem. 

Before stating the Nitsche and Schatz Theorem, let us introduce some known properties of the finite 
element spaces Vjf. 

Assumption 2. Given two fixed concentric spheres Bq and B with Bo cc B c:c 12, there exists 
an ho such that for all 0 < h ^ ho, we have for some R ^ 1 and M > 1; 

B1 For any 0 ^ s < R and s ^ ^ M, for each u e H e (B), there exists rj e Vfi such that 

\\u — v\\s,B < Cl/ 

Moreover, if u e Hq(Bo) then rj can be chosen to satisfy p e Hq(B). 

B2 Let p e ^ 0 °°(Bo) and Uh e Vjf, then there exists p e V p f) Hq(B) such that 

||<puh - p\\i,B ^ C(ip, B, B 0 )h\\uh\\i,B- 

B3 For each h < ho there exists a domain Bh with Bo ccz Bh ccr B such that ifO^s^l^R 
then for all Uh e Vjf we have 

\\uh\\e,B h < Ch s e \\uh\\ s ,B h - 

We now state the following theorem, a key tool in the forthcoming proof of Theorem 1. 


4 





Theorem (Nitsche and Schatz [22]). Let Oo ^1 ^ and let V£ sa ^ s fy Assumption 2. Let 

u £ H\Fl\), let Uh £ Vjf and let q be a nonnegative integer, arbitrary but fixed. Let us suppose that 
u — Uh satisfies 



V{u-u h ) 

■ = 0, 

\/ Vh £ v£ n ffo(fli). 


Jfi 



Then there exists h\ such that if h ^ hi 

we have 


(i) for s = 0. 

, 1 and 1 =5 t ^ M, 




11^ ^/i||s,Qo 

^ C (h l ~ s | 

k tfil + U ~~ UhW-q,^ 

(ii) for 2 ^ s 

^ ^ M and s ^ k < R, 




II^ Uh\\ s ,n 0 ^ C (h ||^||^,ni T h \\u n/j||_ qr,QiJ • 

In this paper, we will actually need a more general version of the assumptions on the approxi¬ 
mation space V£: 

Assumption 3. Given B c II, let be p' ^ 2, there exists an ho such that for all 0 < h ^ ho, we 
have for some R ^ 1 and M > 1: 

B1 For any 0 ^ s ^ R and s ^ ^ M, for each u e H e (B), there exists 77 £ Vjf such that, for 

any finite element T a B, 

\u-v\ s y, T ^Ch d A/p s \u\e,2,T- 

B3 For 0 ^ s ^ ^ R, for all Uh £ Vjf, for any finite element T in the family Th, we have 

\\ u h\\e,p',T Ch d ^ p l l 2S> h s £ ||u/i|| S) 2,t- 

Assumptions B1 and B3 are generalisations of assumptions B1 and B3. They are quite standard 
and satisfied by a wide variety of approximation spaces, including all finite element spaces defined 
on quasi-uniform meshes [ ]. The parameters R and M play respectively the role of the regularity 

and order of approximation of the approximation space Vjf. For example, in the case of Pi-finite 
elements, we have R = 3/2 — e and M = 2. Assumption B2 is less common but also satisfied by a 
wide class of approximation spaces. Actually, for Lagrange and Hermite finite elements, a stronger 
property than assumption B2 is shown in [5]: let 0 ^ s ^ l ^ k, ip £ c tafi(B) and Uh £ Vjf, then 
there exists rj e Vjf such that 

|| Fu h - r]\\ S) B < C(ip)h e ~ s+1 \\u h \\^ B . (5) 

Applied for s = £ = 1, inequality (5) gives assumption B2. 


2.4 Statement of our main results. 

Our main results are Theorems 1 and Theorem 2. The rest of the paper is mostly concerned by the 
proof and the illustration of the following theorems. 
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Theorem 1. Let LIq cc; Qi czcz 0 satisfy Assumption 1, 1 ^ s ^ k. Let u$ be the solution of 
problem (P$) and u k its Galerkin projection onto V k , satisfying (2). Under Assumptions 2 and 3, 
there exists hi such that if 0 < h ^ hi, we have, 

IN - «tf||i,no < C(Q 0 , Lli,Q)h k ^/\ ln/i|. (6) 

In addition, for s ^ 2, if the finite elements are supposed H k -conforming, we have 

IN - uJlU.no < C7(fi 0 ,fii, n)h k+1 - s V\^h\- (7) 

We can show a stronger result than Theorem 1 for Pi-hnite elements. With the same notations 
and under the same assumptions, 

Theorem 2. The Pi-finite element method converges with the order 1 for H 1 (LIq)- norm. More 
precisely: 

\\u$ Ufi ||i,n 0 ^ C (Ho; > Ll)h. 


3 Proof of Theorem 1. 

This section is devoted to the proof of Theorem 1. We first show a weak version of the Aubin- 
Nitsche duality lemma (Lemma 1) and establish a discrete inf-sup condition (Lemma 2). Then, we 
use these results to prove Theorem 1 


3.1 Aubin-Nitsche duality lemma with a singular right-hand side. 

The proof of Theorem 1 is based on Nitsche and Schatz Theorem. In order to estimate the quantity 
N — we will first show a weak version of Aubin-Nitsche Lemma, in the case of Poisson 

Problem with a singular right-hand side. 

Lemma 1. Let f e W~ 1,p {Ll) = (Wq’ p (Ll))', 1 < p < 2, and u e Wo’ p (r2) be the unique solution of 

( —A u = f in LI, 

1 u = 0 on dLl. 


Let Uh £ Vjf be the Galerkin projection of u. For finite elements of order k, letting e = u — Uh, we 
have for all 0 < q ^ k — 1, 

IM|_gr,fi < C/r 9+1 h 2(1/p '” 1/2) [e|i, p ,fi. (8) 

Proof. We aim at estimating, for q ^ 0, the H _9 -norm of the error e: 


e \\-qP 


sup 


Sn e< ? 


0e%°°(fi) IMIU.fi 


(9) 


The error e e Wq ’ p satisfies 


| 

Jo 


Ve • Vv h = 0, \/v h 6 V k . 


Let be i f> e < ^Q°(f l) and let w^ £ H q+2 be the solution of 


—A = cj) in 0, 
w < t > = 0 on dLl. 


6 





In dimension 2, by the Sobolev injections established for instance in [6], H q+2 (Q) c; W 1,p ' ( fl) for 
all p' in [2, +oo[. Thus, for any Wh £ V k , 


f efi 

— 

f eAw+ 

— 

f Ve-Vw* 

— 

f Ve • V(w 0 - w h ) 



Jo 


Jo 


Jn 


We have to estimate \w^ — Wh\iy t fi- It holds 

\ w 0 - ^iC',n = Z ^ 




T 

For all 0 ^ q ^ k — 1 and for all element T in 7thanks to Assumption B1 applied for s = 1, 
l = q + 2, there exists Wh £ Vjf suc ^ as 

Iw' 0 - Wh\i, P ',T C/i 2(1/p,_1/2 )^ +1 |to% +2 ,2,T- (10) 

We number the triangles of the mesh {T*, i = 1, • • • , AT} and we set 

a = (a;)* and b = (^)j, where a* = |u/ - and 6, : = |rc^| 9+2 ,2,T i - 

By (10), we have, for all i in [[1, IV]], 


a, ^ ChWP'-^h^bi. 

We recall the norm equivalence in M. N for 0 < r < s, 

||a;||fs ^ ||x|| (r ^ A rl/,r ~ 1 / s ||a;|j^. 

Remark that here N ~ Ch~ 2 . As 2 < p 7 , we have ||&|L p / ^ \\b\\p- Then, we can write 

\w^ - WhWiyp = ||a|j £p / sS h q+l h 2 ( 1/p '~ l/2:) \\b\\ iP i 

^ h q+1 h 2 ^ p '- l ^\\b\\ P 

< h q+1 h*W-m\vt \ q+2l2l n 

< /l , + i /t 2 (Vp'-i/3)|| 0 || gin . 

Finally, using this estimate in (9), we obtain, for q ^ k — 1, 

||6||_ ff , n <C7W +1 fc 2 (W-i/*)| e | liPin . 


a 


Corollary 1. For finite elements of order k, for any e > 0 small enough, 

||u ,5 ||—fc+i,n ^ C h h \u$ Ufi | i,y>,o• (11) 

Proof. We will apply Lemma 1 to estimate ||tt 5 — u k \\- q ^. This can be done with a suitable choice 
for (p,p'). With s > 0 small enough, let us take: 


P = 


2 

1 + e 


and p 


2 

1 — s' 


In inequality (8): 



Finally, for finite elements of order k, 

II u s - ug II—fc+i,n < Ch k h~ £ \u s - 


= —e. 


HS I l,p,f2■ 


( 12 ) 

(13) 


□ 
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3.2 Estimate of | u$ — u$ 

It remains to estimate the quantity | u$ — w^|i lP ,n by bounding |u^|i iPi n in terms of |r/-<y|(equality 
(15)). To achieve this, we will need the following discrete inf-sup condition. 

Lemma 2. For p and p' defined in (12), we have the discrete inf-sup condition 


mi sup -—j--— - 

UheVf v h eV£ ll u /i||l,p,nF/i||iy. 


^ Ch £ . 


Proof. For e sufficiently small and p and p' defined in (12), the continuous inf-sup condition 


L. Vu • Vu 
inf sup —— > fi > 0 

ueW o P veWfi p ' ll U llbplTlll,p' 

holds for (3 independant of p and p'. It is a consequence of the duality of the two spaces Wq ,p (fl) 
and Wq' v (f2), see [18]. For v e wl' p (Q), let II h,v denote the Ffg-Galerkin projection of v onto Vjf. 
This is well defined since W 1,p (O) c; H ( \(Q). We apply Assumption B3 to IIfor l = s = 1, and 
get 

|iM| iy , n <s ch-w-W) in^llM,!! < Ch-w- w> WlA „ < Cft -2(i/2-i/ P ')n„| iy !! . 

Moreover, for any Uh £ Vjf a W 1,p ( fl), 


\\u h \\i,p,n < C sup 


Sn Vu/i-Vu 


= C sup 
vewfi pl »E^ n iy 


Vttfe • Vllfti) 


II 


< CTT^-W) sup L - vn ^ 

vewy 

^Ch-^-W) sup ^7 Uh ' VVh . 




■ u 7 Ji. 


Finally, thanks to Poincare inequality, and to inequality (13), 

Vu/, • Vvh, 

inf sup -y - - n — y— ^ Ch . 

u h eV h v h evf ll u /il|i,p,n|p/i||i ) p' ) n 


□ 


Then, we can estimate \u$ — u!f |i, p ,n : 
Lemma 3. With p and p' defined in (12), 


| ^5 


< c 



(14) 


Proof. According to Lemma 2, it exists Vh £ Vjf, with ||u/,|| i. p /.n = 1, such that 

h 2{i/2-W)\\ u h\\ lpn ^ c f Vu h . Vvh = c [ Vus . Vvh ^ dlusWivs 

Jn Jn 


I u s 


— ?4li, P ,n < |«5|i, p ,n + |u^|i >P) o < Ch 2 ^ 1//2 1//p ||^||i,p,o- 


So we have 


(15) 









All that remains is to substitute ||uj||i,p,n for the expression established in (4). For p defined as in 

( 12 ), 


ug 






C C_ 
V 2 -p y/e' 


Finally, with (13) and (15), we get 


\us ^5|l,p,Q ^ C . 


Q 


3.3 Proof of Theorem 1. 

We can now prove Theorem 1. 

Proof. The function us is analytic on fii, therefore the quantity is bounded. If we suppose 

s = 1, Nitsche and Schatz Theorem gives, for i = k + 1 and q = k — 1, 

lk<5 - w^lliho ^ C (h k + ||us - • 


Thanks to (11) and (14), 


therefore, taking e = | lnh| 1 , 


ha - us\\-k+i,n ^ Ch 1 


.h~ 2£ 


us - us\\-k+i,n < Ch k y/\ In h\. 


Finally, we get the result of Theorem 1 for s = 1 (inequality (6)): 

ha - ||i,q 0 < Ch k y/\ In h\. 


(16) 


Now, let us fix 2 ^ s ^ fc, Nitsche and Schatz Theorem gives, for £ = k + 1 and q = k — 1, 

ha - ^alkfio < c ( h k+1 ~ s + h 1_s ha - «all-fe+i,ni) • 

So, thanks to (16), we get the second result of Theorem 1 (inequality (7)), 

ha - Ug\\ s ,n 0 < Ch k+1 ~ s yJ\ In h\. 


which ends the proof of Theorem 1. 


□ 


4 Proof of Theorem 2. 

To prove this theorem, we first regularize the right-hand side and prove that in our case the solution 
ug of (Ps) and the solution of the regularized problem are the same on the complementary of a 
neighbourhood of the singularity (Theorem 3). The proof of Theorem 2 is based once again on the 
Nitsche and Schatz Theorem and the observation that the discrete right-hand sides of problem (Ps) 
and the regularized problem are exactely the same, so that the numerical solutions are the same 
too (Lemma 5). 
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4.1 Direct problem and regularized problem. 

Let e > 0, and f e be defined on D by 


fe 


d 

cr(S d _i)e d 


1 B e 


(17) 


where B e = Bfx o,e) and <r(S<i-i) is the Lebesgue measure of the unit sphere in dimension d. The 
parameter e is supposed to be small enough so that B e cc Q. The function f £ is a regularisation 
of de the Dirac distribution Let us consider the following problem: 


(Pe) 



fe in 

0 on dD. 


Since f e £ L 2 (D), it is possible to show that problem ( P £ ) has a unique variational solution u £ in 
iLg(D) f| 5 2 (Q) [16]. We will show the following result: 

Theorem 3. The solution us of (Ps) and the solution u £ of (P £ ) coincide on the closed Ll = Ll\B £ , 
ie, 

ug\„ = u £ I . 
in lo 

The proof is based on the following lemma. 

Lemma 4. Let d £ N*, e > 0, x £ M d , v a function defined on harmonic on B(x,e), and 
f £ L 1 (M rf ) such that 


• f is radial and positive, 

• supp(f) c 5(0,0, e > 0; 

• j f(x) dx = 1. 

Then, f * v(x) = f{y)v(x - y) dy = v(x). 

J R d 

Proof. As supp(f) a 5(0, e), using spherical coordinates, we have: 

f * v(x) = f(r)v(x — ruj)r d ~ 1 du)dr= r d l f{r)[ v(x — ru)duj ) dr. 

Jo J§ d -! Jo \J§i-i J 

Besides, v is harmonic on 5(x,e), so that the mean value property gives, for 0 < r ^ e, 

„d -1 


v(x) = 


a{dB{. 


\x,r)) I 


v (y) d y = tant \\ i 
8B(x,r) a(dB(x,r)) J S d- 


v(x — no) dw, 

Js d -! 


thus 


f * v(x) = f f(r)v(x)a(dB(x,r))dr = v(x) f f f{r)r d 1 do;d?’ = u(x) f f(y) dy = v(x). 

Jo Jo Js ^ 1 Jfi(o,e) 


(S 


Now, let us prove Theorem 3. 
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Proof. First, let us leave out boundary conditions and consider the following problem 

- Au= f e in @'(R d ). (18) 

As —AG = So in 0'(R d ), we can build a function u satisfying (18) as: 

u{x) = fe * G(x) = f f e {y)G(x-y)dy = f f £ (x 0 + y)G(x-x 0 -y) dy = (f £ {x 0 +■) *g)(x- x 0 ). 

jR d v ' 


Moreover, for all x e fi\B £ , G is harmonic on B(x — xo,s), and f e (- + xo) satisfies the assumptions 
of Lemma 4, so that u(x) = (^f £ (xo + -)*G^j(x — xo) = G(x — x o). We conclude that u and G(- —xo) 
have the same trace on dQ, and so u + u, where u is the solution of the Poisson problem (3), is a 
solution of the problem (P e ). By the uniqueness of the solution, we have u £ = u + oj. Finally, for 
all x e fi\B £ , u £ (x) = ug(x). These functions are continuous on Q = £l\B s , this equality is true on 
the closure of ft, which ends the proof of Theorem 3. □ 

Remark 1. Theorem 3 holds for any radial positive function f e L 1 (M (i ) fiL 2 (M <i ) such that 

supp(f) cr B(0,e) and f{x)dx = 1, 

jR d 


taking f £ = /(■ — xq). It is a direct consequence of Lemma f. 

Remark 2. The result is true in dimension 1, taking f £ = —1 i e , where I £ = [xo — e, xq + e] c 
]a,b[= I. In this case, we can easily write down the solutions us and u £ explicitly, 


r 

u$(x) = < 


b — xo 6 — xo 

7 x a 

b — a b — a 

m 

(p 1 

o 

xo — a xo — a 

, x + b 

b — a b — a 

if xe [x 0 ,6]. 

b — xo b — xo 

, x a 

b — a b — a 

if x e [a, xq — 


U e (x) = < 


.T 2 ( xo a + b — 2xo \ 

_ 4i + (2i + 2(6-a) ) X 
a(x o — b) + b(xo — a) Xq + e 2 
2(6 — a) 4s 


if x e [x 0 - s,x 0 + e], 


xo ~ a 
b — a 


x + 


xo - a 

h ~k - 

b — a 


if x e [xq + e, 6]. 


and observe, as shown in Figure 2, that us and u £ coincide outside I £ . 
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Figure 2: Illusration of Theorem 3 in ID. 


4.2 Discretizations of the right-hand sides. 

At this point, we introduce a technical assumption on B £ and the mesh. 

Assumption 4. The domain of definition B e of the function f £ is supposed to satisfy 

Be <= T 0 , 

where Tq denotes the triangle of the mesh which contains the point xq (Figure 3). 



Figure 3: Assumption on B e . 


Remark 3. The parameter e will be chosen to be h/ 10, so it remains to fix a “good” triangle To and 
to build the mesh accordingly, so that Assumption 4 is satisfied. Remark that it is always possible 
to locally modify any given mesh so that it satisfies this assumption. 

Lemma 5. Under Assumption 4, 

4 - 4 . 

where and u^ are respectively the numerical solutions of problems ( P $) and ( P £ ). 
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Proof. Let us write down explicitly the discretized right-hand side F'f associated to the function 
f £ : for all node i and associated test function Vi £ V*, 

i F e)i= I -7^rlB e (a:)ui(x)dx = f —Vi(x) dar, 

Jn<r{ B e) Jb £ cT 0 <r{ B e) 

and Vi is affine (and so harmonic) on Tq, therefore 



if i is a node of the triangle To, 
otherwise. 


We note that Fjf = D h , where D h is the discretized right-hand side vector associated to the Dirac 
mass. That is why, with Ah the Laplacian matrix, 


4-A- 2 [VW 

i node 



Vi = 0. 


□ 


Remark 4. F £ = D h holds as long as B £ c Tq. 

3 ), but Ff A D h , and so ub\ A u £ \ . 

£ ’ £ |n 


Otherwise, we still have us\ n 


= u e 


|n 0 


( Theorem 


4.3 Proof of Theorem 2. 

Theorem 2 can be now proved. 

Proof. First, by triangular inequality, we can write, for s e {0,1}: 

I U S ^ 5 l|s,flo ^ II ^<5 ||s,f2o II ||s,f2o ~b || U$ || . 

Besides, thanks to Theorem 3, we have 


\\ u 5 || s,Oo 0) (19) 

and thanks to Lemma 5, we have 

Il'Utf — u e || s> a 0 = 0- 

Finally we get 

\\u S U$ ||s,£lo ^ ll^e ||s,f2o • (20) 

We will apply Nitsche and Schatz Theorem to e = u £ — v!f. With £ = 2, s = 1, and p = 0, 

ll e l|l,Q 0 ^ C (^ll u e||2,Qi + IMIo.Oi) • (21) 

The domain Pi is a smooth and f £ e L 2 (P), so u £ £ FI 2 {Pi) f| Hq{PI), and then, thanks to inequality 

( 1 ), 

IM|o,Oi < ll e !lo,o < Ch 2 \\u £ \\ 2 ,n. < C7i 2 ||/ e ||o,n. 

As ||/ £ ||o,o can be calculated, 

ll/ ‘ lkn = (l(i lB - (! ' ) ) d!/ 

for e ~ h/ 10 (in order to respect the assumption on B £ ), we get 

||e||o,ni < Ch. (22) 
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Finally, according to Theorem 3, u§\ = u e \ , therefore combining (21) and (22), we get 

|Qi 

II U e - || 1^0 = ||e||i,n 0 < Ch. 

At last, we obtain from inequalities (20) and (23) the expected error estimate, that is 


us - walli,n 0 < Ch. 


( 23 ) 


□ 


Remark 5. In the particular case of the Lagrange finite elements, Koppl and Wohlmuth [19] have 
shown for the lowest order case quasi-optimality and for higher order optimal a priori estimates, 
in the ifi-norm of a subdomain which does not contain xq. The proof is based on Wahlbin-type 
arguments, which are similar to Nitsche and Schatz Theorem (see [27], [28]), and different arguments 
from the ones presented in this paper, like the use of an operator of Scott and Zhang type (see [26]). 
More precisely, if we introduce a domain fl 2 such that LIq czcz fix ccz Q 2 cc LI, xq, f fl 2 , und 
satisfying Assumption 1, with our notations and under the assumptions of Theorem 1, their main 
result is 

II rto o o'! / h 2 \ln(h)\ if k = 1, 

F<5 - ^<5llo,Oi ^ C(i2i,i2 2 ,S2) | hk+1 ifk> 1 24 > 

Using this result, we can easily prove the inequality (6): we apply the Nitsche and Schatz Theorem 
on fio and fii for l = k + 1 and q = 0, 

||w<5 — Us III,o 0 < C (h k + ||«5 — mJ'IIo,^! J • 

Combining it with (24), we finally get 

H - u h 5 ||i,n 0 < C (h k + h k+1 1 ln(h)|) ^ Ch k . 

This result is slightly stronger than inequality (6), but it is limited to the Lagrange finite elements 
and the H l -norm. 


5 Numerical illustrations. 

In this section, we illustrate our theorical results by numerical examples. 

Concentration of the error around the singularity. First, we present one of the computations 
which drew our attention to the fact that the convergence could be better far from the singularity. 
For this example, we define as the unit disk, 

LI = {x = (xi,x 2 ) e M 2 : ||x|| 2 < 1}, 


flo as the portion of 0 

Oo = {x = (xi,x 2 ) e M 2 : 0.2 < ||x|| 2 < 1}, 

and finally xq = (0, 0) the origin. In this case, the exact solution us of problem (P$) is given by 


u(x) = _ ^“ log + 
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When problem (P$) is solved by the Pi-finite element method, the numerical solution v!/ con¬ 
verges to the exact solution ug at the order 1 on the entire domain 17 for the L 2 -norm (see [25]). 
The previous example has shown that the convergence far from the singularity is faster, since the 
order of convergence in this case is 2 (see [19]). The difference of convergence rates for L 2 -norm 
on 17 and 17o let us suppose that the preponderant part of the error is concentrated around the 
singularity, as can be seen in Figures 4, 5, 6, and 7. Indeed, they respectively show the repartition 
of the error for l/h ~ 10,15, 20 and 30. 

Error Error 



Figure 4: Error for l/h ~ 10. 


Figure 5: Error for l/h ^ 15. 



Error 



Figure 6: Error for l/h ~ 20. 


Figure 7: Error for l/h ~ 30. 


Estimated orders of convergence. Figure 8 shows the estimated order of convergence for the 
i7 1 (f7o)-n°rm for the P^-finite element method, where k = 1,2,3 and 4, in dimension 2. The 
convergence far from the singularity (i.e. excluding a neighborhood of the point xo) is the same as 
in the regular case: the P^-finite element method converges at the order k on 17o for the P 1 -norm, 
as proved in this paper with a y/\ ln(/i)| multiplier. 
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Figure 8: Estimated order of convergence for Tf^f^-norm f° r the finite element method P^, k = 

1 , 2 , 3 , 4 . 
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